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1. Introduction. 

The known procedures of generating the dynamics of a conservative mechanical system from the 
action functional, the Lagrangian or the Hamiltonian do not apply directly to dissipative systems. In 
preparation for a variational formulation of dynamics with dissipation we are adapting in the present 
note concepts borrowed from convex analysis to differential geometric setting. We are applying these 
concepts to examples of static systems. Holonomic and non holonomic constraints as well as partial 
control are considered. We believe that the principle of virtual work applied to static systems is a 
simple fundamental model of all variational principles of classical physics. 

The paper is organized as follows. In Section 2 we state the principle of virtual work which is 
incorporated in the definition of the constitutive set. In Sections 3-6 we adapt concepts of convex 
analysis to static systems. In Sections 7-8 we consider the case of partially controlled potential systems. 
A generalization to the case of dissipative systems leads to the notion of generating families of forms 
and their critical points (Sections 9-10). 

Propositions and theorems stated in this note without proofs are found in [2]. 

2. The principle of virtual work and examples of static systems. 

The principle of virtual work well known in statics of mechanical systems is a master model for all 
variational principles of classical physics. We state a simple version of the principle. The configuration 
space of a static system is a set Q assumed here to be a differential manifold. The system is represented 
by its internal energy 

U:C°^R (1) 

defined on a constraint set C° C Q assumed to be a submanifold. The principle of virtual work is 
incorporated in the definition 

S = {/ e T*Q ; q = n Q (f) e C°, V, 9eT?c o (dU,8q) = (f,8q)} (2) 

of the constitutive set. The constitutive set characterizes the response of a static system to control by 
external forces. 

The configuration spaces in the examples of static systems in the present section are constructed 
from the Euclidean affine space M of three dimensions modeled on a vector space V. The Euclidean 
metric is represented by a metric tensor g : V — > V* . 
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Example 1. Let the configuration space of a material point be the afHne space Q = M. The 
space SQ = V is the model space of Q and F = V* is the dual of the model space. The point with 
configuration q G Q is tied with a rigid rod of length a to a fixed point with configuration q . The 
configuration q is constrained to the sphere 

C° = {q€Q; \\q-qo\\=a}. (3) 

With the internal energy U = the constitutive set is the set 

S= {(?,/) GQxF; \\q-qo\\=a, f = a~ 2 (f,q - q )g(q - q )} . (4) 

▲ 

A more general static system is characterized by virtual work function 

o~ : C 1 — > R (5) 

defined on a constraint set C 1 C TQ. For each geC° = tq(C 1 ) the set C\ = C 1 (~l T g Q is a cone in 
the sense that if Sq G C*, then AJg G C* for each A ^ 0. The virtual work function is homogeneous in 
the sense that a(XSq) = Xa(6q) if A ^ 0. The inclusion C 1 C TC° is usually verified. Constraints are 
said to be holonomic if C 1 = TC°. The tangent set TC of any set C C Q is defined by 

TC = {Sq e TQ ; there is a curve 7 : R — > Q such that t7(0) = <5g and 7(3) G C if s ^ 0} (6) 

Holonomic constraints are fully represented by the constraint set C . Constraints are said to be bilateral 
at q e C° if (5g e implies that — Sq eCj. If constraints are not bilateral at q G C°, then they are 
said to be unilateral. Constraints in the above example are holonomic and bilateral. 
A more general version of the principle of virtual work is incorporated in the definition 

S = {/ G T*Q ; q = n Q (f) G C°, V 5(?eC , <r(5q) > (f,Sq)} (7) 

of the constitutive set. 

Example 2. Let the configuration space of a material point be the affine space Q = M. We use 
definitions and symbols introduced in Example 1. The material point is not constrained and is subject 
to isotropic static friction. The virtual work is the function 

o : Q x SQ -> R : (q, Sq) » p(q)\\Sq\\ = p(q)y/ '(g(8q),6q) . (8) 

The set 

S = {(«,/) G Q x F ; y 5qeSQ P (q)\\5q\\ > (f,Sq)} (9) 

is the constitutive set. 

Let (q, f) G S. By setting Sq = g~ 1 (f) in the inequality 

p(q)\\SQ\\>{f,Sq) (10) 

we obtain the inequality 

P(q)\\f\\ > ll/ll 2 - (ii) 

Hence 

SC {(«,/) GQxF; ll/H <p(g)}. (12) 

Let (q, f) satisfy the inequality 

ll/ll (13) 

T^ 1 tic rol^it ion 

(/,^)<|(/,^)K||/||||^||<p(g)||^|| (14) 

is derived from the Schwarz inequality 

\{fM)\^\\f\\¥q\V (15) 

We have shown that 

S = {(q,f)GQxF- \\f\\^p(q)}. (16) 

▲ 
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Passages similar to that from (9) to (16) will appear in several examples. 

Example 3. The configuration space of a skate is the set Q = M x D, where D is the projective 
space of directions in M. We use the Euclidean metric in M to identify the space D with the unit 
circle 

fl = {i?e^;(#),t?)=l}. (17) 
Virtual displacements are elements of the space TQ = M x V x TD, where 

TD = {(<?, 50) g D x V ; J0) = 0} . (18) 

The skate is a system with non holonomic constraints. The set C° is the entire space Q. The constraint 
consists in restricting virtual displacements in M to those parallel to the direction specified by an 
element of D. Thus 

C 1 = {(a;, Sx, 0, <W) e ; 3 AeR = Atf} . (19) 
The constitutive set is a subset of the space T*Q = M x V* xT*D, where 

T*D = {(i?,t) e L> x ; (r,i?) = 0} . (20) 

is the space chosen as the dual of TD. 

Let the skate be subject to friction represented by a non negative function p : Q -»R, The virtual 
work is the function 

a-.C 1 -^R:(x,5x,$,5ti) ^ p(x,d)\\5x\\ = p(x, ti)y/(g(8x), 6x). (21) 

The set 

S={(i,/,tf,r)gT*Q;V WliW ^MHNI > </,fe) + (r,W)} (22) 
is the constitutive set. The equality r = is obtained by setting <5x = in the inequality 

p{x,0)\\5x\\>{f,5x) + (T,50) (23) 

with arbitrary By setting <5x = Ai9 we arrive at the inequality 

p(x,&)\\\>\(f,0) (24) 

for each A 6 I. The inequality must be satisfied for A = (/,#). Hence p(;r, $)|(/, > {f,$) 2 and 
|(/,<)|<P(M). K K/,0>|<p(aM?),then 

p(z,0)|A| > |A||</,0)| ^ (/, Ai?) (25) 

for each Ael. Hence, the virtual work principle is satisfied. 
In conclusion we obtain the expression 

S={(xJ,$,r)eT*Q: |</,0)| < p(x,#),r = 0} (26) 

for the constitutive set of the system. ▲ 

Example 4. The present example gives a formal description of experiments performed by Coulomb 
in his study of static friction. Notational conventions of Example 1 will be used. Let a material point 
with configuration q £ Q = M be constrained to the set 

C° = {qeQ; (g(k),q-q ) ^0}, (27) 
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where qo is a point in Q and fc G V is a unit vector. The boundary 

dC° = {geQ; ( 3 (fc), g - go } = 0} (28) 

is a plane passing through q and orthogonal to fc. In its displacements along the boundary the point 
encounters friction proportional to the component of the external force pressing the point against 
the boundary. The system is characterized by the virtual work function a = defined on the non 
holonomic constraint 

C 1 = {(q, 5q)eQxSQ; (g(k),q - q ) > 0, 

(g(k), Sq) > py/\\6q\\*-{g(k),6q)* if (g(k),q - q ) = o} , (29) 

where p > is the coefficient of friction. 

The principle of virtual work defines the constitutive set 

S = {(q, /) G Q x F ■ y 5qee5Q if (q, Sq) G C\ then (/, Sq) ^ o} . (30) 

If the material point is not on the boundary, then (g(k), q — qo) > 0. The virtual displacements are 
not restricted and (q, f) G S if and only if / = 0. 

If the material point is on the boundary, then (g(k),q — qo) = 0. We show that in this case (q, /) G S 
if and only if the inequality 

V\\f\\ 2 -(f,k) 2 +p(f,k)^0 (31) 

is satisfied. The inequality ||/|| 2 — (/, fc) 2 ^ is guaranteed by the Schwarz inequality | (/, k) | ^ ||/||||fc||- 
Let (q,f) G S and let ||/|| 2 - (/, fc} 2 = 0. This will be the case when / = -\\f\\g(k). The case 
/ = ll/llffC 1 ) with ll/H ^ is excluded by the virtual work principle. The inequality (31) is satisfied 
since 

y/\\f\\ 2 -(f,k) 2 + p(f,k) = -p\\f\\^0. (32) 
We analyse the case of (q, f) G S and ||/|| 2 — (/, fc) 2 > 0. The virtual displacement (q, Sq) with 



Sq = 9~\f) ~ (/, k)k + pVll/ll 2 - if, k) 2 k (33) 

is in C 1 since 

(g(k),Sq)=p^\\fP-(f,k)*. (34) 
From the principle of virtual work and from 



</, Sq) = ll/f - (/, fc) 2 + PvW~W(/, k) (35) 

it follows that 

ll/f - (f,k) 2 + p^f¥^JKW(f,k) <0. (36) 

The inequality (31) is obtained since ||/|| 2 - (/, fc) 2 > 0. 

Let the point be on the boundary and let the inequality (31) be satisfied. The Schwarz inequality 



\(g(u),v) - (g(k),u){g(k),v)\ < V\H\ 2 W), «) 2 V\M 2 - (g(k), «) 2 (37) 
for the bilinear symmetric form 

(u,v) i ► 0/(«),«> - (g(k),u)(g(k),v) (38) 
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applied to the pair (g 1 (f),Sq) leads to the inequality 

(/, Sq) - (/, k)(g(k),6q) < V\\.f\\ 2 -(f,k)W\\Sq\\ 2 -(g(k),5q)2. (39) 
If y/\\f\\ 2 -(f,k)* + p(f, k) < 0, and (g(k), Sq) > p^/WW^WhW, then 

VU¥~lfWV¥q¥~j9WW ^ -P(f, k)V\\Sq\\ 2 -(g(k),6q)Z ^ -(/, k)(g(k), Sq). (40) 

It follows that (f,Sq) ^ 0. Hence, (q, f) is in the constitutive set S. 
We have shown that 

S = {(q, f) G Q X F ; < 5 (fc), g - go) ^ 0, / = if q - q Q ) > 

and \/||/|| 2 - </, fc) 2 + P</, fc> < if ( 9 (k), q-q o ) = 0}. (41) 

We will interpret the friction coefficient p in terms of the components 8q x = (g(k),5q)k and Sq" = 
Sq — (g(k), 5q)k orthogonal and parallel to the surface dC°. The inequality 

(g(k), Sq) > pvfsFWW (42) 

translates into the inequality 

\\Sq-\\^p\\Sq"\\ (43) 

for the norms 

\\Sq^\\ = (g(k),Sq) and \\5q" || = ^INI 2 - (g(k), 8q)\ (44) 

Let q be on the boundary of C° and let d be the angle between a displacement and fc. The pair 
(q, Sq) is in C 1 if i? = or cot ^ p. A 

3. Convex sets and functions. 

Let V be a vector space of finite dimension m. The convex hull of a subset A C V is the set 

cnv A = { v e V ' 3 01;02; ... )0 , eA , Al ,A 2l ... ,A fc e[o,i] Sti = 1 and v = Eti ■ (45) 

The inclusion A c cnv A is always satisfied. The set A is said to be convex if A = cnv A. The convex 
hull cnv A is the smallest convex subset of V containing A. The affine hull of A is the set 

aff A = { v e V : 3 01;02; ... , afceA , Al , A2i ... , Afc£K Eti = 1 and « = Eti A *°*} • (46) 

It is the smallest affine subspace of V containing A. 

Let A C V be a convex subset. The set aff A is a topological space. The relative interior rintA of 
A is the topological interior of A in aff A. It is the biggest subset of A open in aff A. 

Proposition 1. The relative interior rint A of a non empty convex set A C V is not empty. 

Let o~ : C — > K be a function defined on a subset C C V. The overgraph of cr is the set 

ovgr a = {{v,r) eV xR; v eC,r ^ er(v)} . (47) 

The function a is said to be convex if ovgr a is a convex subset of V x M. The domain C of a convex 
function is convex. The function a is said to be concave if the function — a is convex. A convex function 
a is said to be closed if ovgr a is closed in V x K. A concave function cr is said to be closed if the 
function —a is closed. 
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Given a static system characterized by constraint sets C° and C 1 and a virtual work function 
a : C 1 ->lwe introduce sets 

Cl = {8q€5Q; (q,5q) € C 1 } (48) 

for each geC" and functions 

a q :C\ Sq » a(q, Sq). (49) 

The static system is considered convex if the functions a q are convex. The three systems in the typical 
examples are convex. We list the objects S q and a q in the three examples. 

(1) Example 1: 

C\ = {Sq € SQ ; (g(q - q ),Sq) = 0} , (50) 
a q = 0. (51) 



(2) Example 2: 



(3) Example 3: 



C\ = SQ, (52) 
a q :SQ^:Sq^p(q)\\Sq\\. (53) 



C(x,tf) = { ( fa > Stf)eVx T$D ; 3 AeK ^ - At?} , (54) 
(4) Example 4: 

= {Sq e SQ ; (g(k),Sq) > Py /\\6q\\* - (g(k),6q)) 2 ^ (flW, 9 - 5b> = o} , (56) 

<Tg = 0. (57) 



4. Separation of sets. 

A non zero affine function h : V — > R is said to separate non empty subsets A C V and _B C 1^ if 
> and h\B ^ 0. The separation is said to be strong if h\A > and h\B < 0. If a function h 
separates sets A and B, the level set 7? = /i _1 (0)) of the function is said to separate the sets A and 
£>. If the function h separates the sets A and B strongly, then the level set H = /i _1 (0)) is said to 
separate the sets A and B strongly. 

The following theorem is known as the Separation Theorem. 

Theorem 1. Let A and B be convex subsets ofV. If 

rint A n rint .8 = 0, (58) 
then there is an afhnc function h separating the sets A and B. 

Theorem 2. Let A c V and B c V be convex, closed and non empty. If A is compact and AnB = 0, 
then there is an affine function h strongly separating the sets A and B. 

5. The Legendre transformation in a vector space, the homogeneous case. 

Let 

cr : C -» R (59) 

be a function defined on a subset C of a vector space V. The set C is assumed to be a cone in the sense 
that if v € C, then kv e C for each fc ^ 0. The function cr is assumed to be positive homogeneous in 
the sense that cr(fcw) = ka(v) for each fc ^ 0. The set 

S={f€V*;V vf . c a{v)-{f,v)>0} (60) 

is called the Legendre transform of a. 
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Proposition 2. The Legendre transform S of a positive homogeneous function a : C —> R is convex 
and closed. 

PROOF: If /i e S and / 2 E 5, then for each v E C and each s 6 [0, 1] we have 

a(v) - <(1 - + sfc,v) = (1 - s)a(v) - <(1 - s)f u v) + sa(v) - (sf 2 ,v) 

= (1 - s)(a(v) - (/!,«» + s(tr(v) - (/ 2 , «)) ^ (61) 

since er(v) — (fi,v) ^ and a{v) — (f 2 , v) ^ 0. It follows that (1 — s)fi + sf 2 E S. Hence, S is convex. 
Continuity of the mapping ( • , v): V* — > M: / h- > (/, v) for each v E V implies that S is closed. ■ 

The constitutive sets S derived from the virtual work functions a in the examples of static systems 
are obtained by applying of the Legendre transformation to functions a q defined in (49). The Legendre 
transforms S q are then combined. The constitutive sets are the unions 

U M x S *- ( 62 ) 
qeQ 

We list the Legendre transforms obtained in the three examples. 

(1) Example 1: 

S q = {fEF; f = a- 2 (f,q-q Q )g(q-q )}. (63) 

(2) Example 2: 

S q = {fEF;\\f\\^p(q)}. (64) 

(3) Example 3: 

S {x ,#) = {(f,r) EV*x T*D ; ^ p(x),T = 0} . (65) 

(4) Example 4: 

S q = {fEF; V\\f\\ 2 -(f,k) 2 + p(f, k) < if (g(k), q - q ) = o} , (66) 

6. The inverse Legendre transformation. 

Let S be a subset of V* . We introduce the set 

C = {vEV; sup fes (f,v) <oo} (67) 

and the function 

(7 : C — > M: ?) i — > supygg(/, v). (68) 
Note that the set C is a cone and a is positive homogeneous. 
PROPOSITION 3. The set C is convex. The function a is convex and closed. 
PROOF: For v\ E C, v 2 E C, and s E [0, 1] the combination svi + (1 — s)v 2 is in C since 

sup /eS (/,st>i + (1 - s)v 2 ) <sup /eS s(/,Di) +sup /es (1 - s)(f,v 2 ) 

= ssu P/es(./>i) + (! - s ) su P/es(./>2) < oo. (69) 

It follows that C is convex. Moreover, 

<t(sui + (1 - s)v 2 ) scr(ui) + (1 - s)ct(w 2 )- (70) 

Hence, a is convex. 
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Let (v,r) be the limit of a sequence (v n ,r n ) of elements of ovgrcr. The relation 

r n > <j(v n ) = sup feS (f, v n ) (71) 

is a consequence of the construction of o and the definition of the overgraph. It follows from this 
relation that for each e > the inequalities 

(f,v n ) ^r n <r + s (72) 

hold for all / e 5 and sufficiently large n. Convergence of (f,v n ) to (/, v) implies that (/, v) ^ r + e 
for each / G S and each e > 0. Hence, (/, v) ^ r and v e C. Moreover 

su P/es(./» < r + £ ( 73 ) 

for each / G S and each e > 0. Hence, <j(v) < r and (w, r) e ovgrtr. We have shown that o is closed 
since ovgrcr is closed. ■ 

Rockafellar [2] calls the function u the support function for the set S. The following theorem is based 
on the results established in [2]. 

Theorem 3. The Legendre transformation and the inverse Legendre transformation establish a one 
to one correspondence between positive homogeneous closed convex functions defined on cones in V 
and non empty closed convex subsets ofV*. 

It follows from Theorem 3 that the constitutive set provides a complete characterization of a convex 
static system. Systems presented in the three examples above are convex. Each of the systems is 
equally well characterized by the virtual work function a or by the constitutive set S. 

7. Generating families and the principle of virtual work for partially controlled systems. 

There are two configuration spaces involved in the description of a partially controlled system. 
There is the internal configuration space Q and the control configuration space Q. The two spaces are 
components of a differential fibration 



Q 



Q 



(74) 



There is the internal energy function 

U:Q^R (75) 

interpreted as a family of functions defined on fibres of the fibration n. The symbol (U, rf) is used to 
denote this family. 

A generating family (U, n) generates the constitutive set 



S = j/ G T*Q ; 3_ e ^ m = n Q (f) y s - eJ _Q (dU,Sq) = (/, T V (Sq))} 



(76) 



of a partially controlled system. 
We denote by \IQ the subbundlc 



{Sq G TQ ; Tv(Sq) - 0} (77) 



of vertical vectors. The set 

Cr(U, ri) = {qeQ; (dU, Sq) = for each Sq £ V ? Q} (78) 

is called the critical set of the family. If q satisfies the conditions stated in the definition of 5, then 
the equality (dU(q),5q) = is obtained with Sq = and any vertical vector Sq £ Q-^. It follows that 
q£Cr(U,r ] ). 

There is a mapping 

n(U, rj) : Cr{U , rf) — ► J*Q (79) 

characterized by 

(K(U,r,)(q),Sq) = (dU,6q) (80) 

for each 5q £ T^^Q and each Sq £ TqQ such that Tr](Sq) = Sq. The constitutive set is the image of 
K{U,r,). 

Example 5. A material point with configuration q 1 in the affine space Q — M is connected to a 
fixed point qo with a rigid rod of length a. The model space V of Q will be denoted by SQ and the 
dual space V* will be denoted by F. The set 

D = {■& £ V ; (5(1?), 1?) = 1} (81) 

will be used as the configuration space of the point q 1 . Its actual configuration of q' in Q is given by 
q' = qo + a-d. A second material point with configuration q is tied elastically to q' with a spring of 
spring constant k. The configuration space is the product Q = Q x D. The function 

U:Q^R:(q,#)~ ||k - (g + ad) f (82) 

is the internal energy. Spaces TD and T*D are defined as in Example 3. Virtual displacements are 
elements of the space TQ = Q x SQ x TD. The constitutive set is a subset of the space T*Q = 
Q x F x J*D. The set 

S={(q, f, tf, t) £ T*Q ; / = kg(q - (q + ad)), r = -ka (g(q - q ) - (g(q - «,),%(<?))} (83) 
is the constitutive set. 

If the configuration -d is not controlled, then r = and we have an example of a partially controlled 
system. The space Q is the internal configuration space and the control configuration space is the 
space Q = M. The canonical projection 

V :Q ^Q: (q,0) » q (84) 
is the relation between the two spaces Q and Q. The internal energy defines a family (U, 77). The set 
Cr(U, V ) = £Q- d £ D if q = q , 1? = ±(q - q Q ) \\q - q \\~ 1 } . (85) 

is the critical set. The description 



K(U, V )(q,ti) 



—kag(d) if q = qo and d £ D 

kg(q - <lo) (1 - a\\q - q-oir 1 ) if q + qo and 1? = +(g - q )|k - QolT 1 (86) 
. M<7 - ?o) (1 + a||<? - fZolT 1 ) if 9 7^90 and 1? = -(q- q )\\q - <j n || _1 



of the mapping k(U, rj) : Cr(U, 77) — > T*Q = Q x is derived from the general expression 

K (U,r 1 )(q,ti) = kg(q-(qo + a#)). (87) 
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The set 



S={(q,f)eT*Q; \\f\\ = ka if q = q , 

f = fe(i ±a||g - 90 II" 1 ) g(q - go) if g^go} 

is the constitutive set. It is the image of k(U, rj). 

Note that the critical set is not the image of a section of rj. For each control configuration q we have 
two different internal equilibrium configurations (q, i?) if q ^ qo and an infinity of internal equilibrium 
configurations if q = qo- The external force necessary to maintain the control configuration q depends 
on the internal configuration. Thus even if the internal configuration is not directly observed its 
presence can not be ignored. ▲ 

8. Reduction of generating families. 

Let (U,rf) be a family generating the set (76). We have the following obvious proposition. 

Proposition 4. Let q e Cr(U, n). The single point set 

% = {/ G T*Q ; n Q (f) = V (q) dU(Sq) = (/, J V (Sq))} ■ (89) 

is represented in the form 

S<i = {/ G T*Q ; ttq(/) = V (q) V 59eT?j( _ )Q ^(<5g) = (/, 5q)} , (90) 

where 

o-? : T,(3)Q -» R: ^ ^ dZ7(^), <5g £ T ? Q, T»/(5g) = 5g. (91) 

It follows from the above proposition that if Cr(U, rf) is the image of a section £ : Q — > Q of the 
fibration then the family (£7, 77) generating the set S in (76) can be replaced by the function 

<T:TQ->R:(6q)~a c{TQ{Sq)) (6q), (92) 

where cr^( TQ (,5 9 )) is the function u-q defined in the the formula (91) with q = C( t q{^q))- It i s obvious 
that a — d(U o (). Thus the set S is generated by the function U = U o (. 

Example 6. Three material points with configurations qo, q, and q' in the affine space M are 
interconnected with springs with spring constants k, k', and k" . The point qo is fixed and not controlled. 
The two points q and q' are not constrained. The configuration q' is not controlled. The internal 
configuration space is the affine space Q = M x M of internal configurations q — (q, q') modelled on 
6Q = V © V. The control configuration space is the space Q — M of controlled configurations q and 
SQ is the model space. The dual of 5Q will be denoted by F. The canonical projection 

n:Q^Q:q=(q 7 q')^q (93) 

is the relation between the two spaces. The internal energy is the function 

U: Q - R:q = (q,q') ~ ^\\q - q \\ 2 + j\\q' - qo\\ 2 + y ||g' - <?|| 2 . (94) 

The internal energy defines a a family (U, rj) of functions on fibres of the projection 77. The critical set 
Cr(U, n) = {q= (q, q') £ Q ; k'g(q' - q ) + k"g(q' - q) = 0} (95) 
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of the family is the image of the section 

C : Q -> Q : q = q h-> (q, go + k"(k' + fc")" 1 ^ - %)) (96) 

of the projection r\. 

The constitutive set is the set 

S = {(q, /) e Q x F ; / = + ^ + - . (97) 

Note that the presence of the material point with configuration q' can be ignored. This is due to the 
fact that the critical set is the image of a section of the projection 77. The constitutive set is generated 
by the reduced internal energy function 

U = U o C : Q -» R: q -» \ ™ \ f ± \\q go f ■ (98) 

▲ 

9. Generating families of forms. 

A generating family consists of a differential fibration 



Q 



v 



(99) 



and a form ct : T Q — ► R. We refer to functions on TQ as forms in order to distinguish such functions 
from functions on Q. We do not imply that the forms are linear. The form <r is differentiable on TQ 
with the image of the zero section removed. It is positive homogeneous and convex on fibres of the 
tangent fibration t-q-.TQ — > Q. We denote by VQ the subbundle 

{SqeTQ; T v {Sq) = 0} (100) 

of vertical vectors. The form <r defines a family (<r, 77) of forms er,j associated with fibres of the fibration 
77. Each form a q is the restriction of the form a to the set 

{SqeTQ; V (TQ(Sq))=q}. (101) 

The set 

Cr(a, 77) = {q G Q ; cr(<5g) > for each <5g G V,Q} (102) 

is called the critical set of the family. 

A generating family (cf, 77) generates the set 



S = {/ G T*Q ; q = ttq(/) G Q, 3_ £ _ if <5<Z € T q Q, 



Sq G T ? Q, and Jn(Sq) = Sq, then cr(<5g) ^ (/, <fy) j . (103) 

If q satisfies the conditions stated in the definition of S, then the inequality a(Sq) ^ is obtained 
with Sq = and any vertical vector Sq G VqQ. It follows that q G Cr(a, 77). Consequently, 

S= (J (104) 

where 

= {/ eT*Q:q = ttq(/) = 77(g), if 5g G T g Q, ^ G T ? Q 

and Tjy(5g) = 5gr, then a(Sq) > (/, <5g}} . (105) 
Proposition 5. Let q e Cr(a, 77). The set Sq is not empty. 
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Proof: It follows from the Separation Theorem that there is a hyperplane in TqQ x R separating the 
subspace VqQ x {0} and the overgraph ovgrcfg of a-q = a\TqQ. This hyperplane is the graph of a linear 
function h on TqQ. Note that a(Sq) > h(5q) and h(Sq') = if Sq' £ V^Q. It follows that the function 

f:T q Q ^R:5q^> h{Sq), Sq = Tr](Sq) (106) 

is well defined, and f £ Sq. ■ 
The relation 

K (a,r)):Cr(a,r)) ^ T*Q (107) 

defined by 

graphK^.T?) = {(?,/) £ Cr(a,r]) x T*Q ; / = Sq} (108) 

generalizes the mapping k(U, rj) introduced in Section 7. We will refer to the set Sq as the contribution 
to the constitutive set S from the critical point q. 

Example 7. A material point with configuration q in the affine space Q = M is tied to a point 
q' £ M with a spring of spring constant k. The point q' is subject to friction and left free. The internal 
configuration space is the affine space Q = M x M of internal configurations q = (q, q') modelled on 
5 Q = V © V. The control configuration space is the space Q = M of controlled configurations q and 
SQ is the model space. The dual of SQ will be denoted by F. The canonical projection 

V :Q^Q:(q,q')^q (109) 

is the relation between the two configuration spaces. The virtual work is the form 

a : Q X 5 Q -> R : (q, q' , Sq, Sq') ^ k{g(q' - q), Sq' - Sq) + p\\Sq'\\ (110) 

This form together with the projection 77 define a generating family for the constitutive set. The set 

Cr(a, v) = {q= (q, <?') £~Q ; V Sq , eV k(g{q' - q), Sq') + p\\Sq'\\ > o} 

= {<j=(q,q')£Q; P>k\\q'-q\\} (111) 

is the critical set. 

The contribution S^q') from the critical point (q, q') to the constitutive set S is the set 

S{q,g>) = {(?./) e Q x F ; Va, e4Q (/,*«) «S g utf v cr(q,q',6q, Sq') 
= {(q,.f)eQxF; \f Sqe5Q (/, tfg) < %(<? - </), 

= {(q,f) xF; / = -«')}, (112) 



since 



due to 



inf <7(<z, <?', Sg, V) = inf (k(g(q' - q), Sq' - Sq) + p\\6q'\\) = -k(g(q' - q), Sq) (113) 

dq'EV dq'EV 



k{g(q' -q),Sq')+p\\Sq'\\ ^ 0. (114) 

The expression 

S= (J S M) ={(q,f)£QxF;\\f\\^p} (115) 

(g,g')£Cr(?,r)) 

is obtained for the constitutive set. Neither the configuration q' nor the spring constant k appear in 
this expression. For each configuration q £ Q the configuration q' is in the set 

{q'£M- {q,q')£Cr(o,f 1 )} = {q l £M ■ k\\q' - q\\ < p} , (116) 
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Example 8. A point with configuration q in Q = M is tied to a point with configuration q' € M 
with a spring of spring constant k. The point q' is in turn tied to a fixed point qo G M with a spring of 
spring constant k! . The point q' is subject to friction and left free. We use the definitions and symbols 
introduced in Example 7. The virtual work function of the system is the form 

a : Q x SQ -> R : (q, q', Sq, 5q') ~ k'(g(q' - q ), Sq') + k(g(q - q'),6q - 5q') + p\\Sq'\\. (117) 

It defines a generating family for the costitutive set together with the projection r\. The set 

Cr(a, T]) = {q= (q, q') EQ ; ^ Sq , eV k'(g(q' - q ), Sq') + k{g(q' - q), Sq') + p\\Sq'\\ > o} 

= {9=(w')eQ; \\k'g(q' -q ) + kg(q' -q)\\^p} (118) 

is the critical set. 

If (q, q') is a critical point, then 



due to 



inf a{q, q', Sq, Sq') = k(g(q - q'), Sq) 

Sq' 



k'{g(q' - qo), Sq') + k{g(q' - q), Sq) + p\\Sq'\\ > 0, 



(119) 



(120) 



and the contribution S^ q<q ^ to the constitutive set consists of a single covector / = kg(q — q'). The 
inequality 

\\k'g(q'-qo) + kg(q'-q)\\^p (121) 



is equivalent to 



It follows that 



is the constitutive set. 



kk! 

k 9(i - 1') - * + - qo) 



kp 
k + k' 



= {(<?, 



f)eQxF; 



kk' 
k + k 



-g(q - qo) 



kp 
k + k' 



(122) 



(123) 



Example 9. We return to the static system of Example 5 modifying the system by subjecting the 
point with configuration q' to isotropic homogeneous static friction. The function 



a:QxSQ -> R : (g, tf, <Sg, i-> k(g(q - (q + a&)),5q) - ka(g(q - (q + ai?)),&?) + pa 



(124) 



together with the projection $ define a generating family of forms. The critical set of the family is the 
set 

Cr(a, V ) = {(q,#)eQ; 0e D, V WeTD - ka{g(q - q ), 50) + pa\\ 80\\ > 0} . 
The inequality 



k\\q -qo- (g(q - qo), = ky/\\q - q \\ 2 - (g(q - q ),#) 2 ^ p 
is derived from the inequality 

-ka(g(q' - q ) - (g(q - q ), 60) + pa\\50\\ = -ka(g(q' - q ),6ti) + pa\\50\\ > 

by setting 



(125) 
(126) 



If 



50 = q-q o -te(q-q o ),0)0. 

k\\q -qo - (g(q - qo),ti)d\\ ^ p, 
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(127) 

(128) 
(129) 



then the inequality (127) follows from the Schwarz inequality 

k(g(q - Qo) - (.9(9 - 9o),t?)3W,W) < k\\q - q - (g(q - %), 0>0|p0||. (130) 
We have shown that 

Cr(a,n) = {(q,ti) €~Q ; tfeD, ^vlFrf^IsFftPP ^} ■ (131) 
The contribution to the constitutive set from (g, 1?) G Cr(a, rf) is the single covector 

/ = *ff(«-(«b-<M?))- (132) 

The set 

S={(?,/)eQxf; 3„ eD W\\q - q \\ 2 - (g(q - <Zo),^> 2 < P , / = fc<?(«/ ~ («0 - <"?))} (133) 
is the constitutive set of the system. 

▲ 

10. Reduction of generating families. 

Proposition 6. Let a : A — > R be a convex function defined on a convex subset A of a vector space 
V. Let W be a subspace ofV, let 

pr:V^V/W:v^[v] (134) 
be the canonical projection, and let B = pr(A). The function 

p : B -> R: [v] ^ inf a(v + w) (135) 

is convex if well defined in the sense that the defining formula (135) assigns to the function finite values. 

Proof: Suppose that p is well defined. Let [v], [v 1 ] s E/V and ^ s ^ 1 For each e > there exist 
vectors v e <E [v],v' e G [v'] such that cr(u e ) ^ p([v]) + £ and cr(^) ^ p([v'}) + e. We have the inequalities 

s + sp([v}) + (1 - s)p([v'}) > *7(v e ) + (1 - S )<7«) 

^ cr(s-y £ + (1 - s)v' £ ) 

>p(s[v] + (l-s)[v'}) (136) 

for each e > 0. It follows that 

sp([v}) + (1 - s)p([v'} > p(s[v] + (1 - s)k/]). (137) 

Hence, p is convex g 
Similarly, for a concave function cr : V — > R, the function 

p : V/W -> R: [u] sup cr(w + 10) (138) 

is concave. 

If the set A is a cone, the set B is a cone as well. If a is positive homogeneous, then p is positive 
homogeneous. If the function a is convex and closed, it may happen that the function p is not closed. 
Let i] : Q — > Q be a differential fibration and let (a 7 , 77) be a generating family. 
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Proposition 7. If q e Cr(a,n), then the set 



Sq = {/ G T*Q ; tt q (/) = 77(g) V 5 _ £T ^ a(ig) > (/, T„(ig)>} (139) 
is represented in the form 

Sq = {/ G T*Q ; ttq(/) = 77(g) V 59eT?)(?)Q ^(£g) > (/, <5g>} , (140) 

where 

Cg : T^(q)Q > R: (5g i > inf o^g), <5g G TgQ, T?7(<5g) = Sq. (141) 

Proof: We show that the function o-q is well defined. As in Proposition 5 we have the graph of a 
linear function h on TqQ separating the subspace V^-Q x {0} and the overgraph ovgrCTq of Wq — a\TgQ. 
It follows from the Separation Theorem that there is a hyperplane in TqQ x K separating the subspace 
^JqQ x {0} and the overgraph ovgrer^ oia-q = a\TgQ. This hyperplane is the graph of a linear function 
h on TqQ. Note that a(Sq) ^ h(Sq) and h(5q') = if Sq' G VqQ. It follows that o-q is well defined since 

inf _a{Sq + Sq) ^ h(5q). (142) 
<5?'eV ¥ Q 

If Sq e Tjj^Q and / G T*^Q, the inequality 

<rq{8q)>(f,6q) (143) 
for each 5g G TqQ such that T?7(<5g) = Sq is equivalent to 

inf a(Sq) > (/, <5g>, 5g G T ¥ Q, Tr?(5g) = <5g. (144) 

Sq 

■ 

It follows from Proposition 7 that if CV(cf, 77) is the image of a section £ : Q — > Q of the fibration 77 
then the family (a, 77) generating the set (103) can be replaced by the generating form 

<t:TQ ->R:(5q) ^ CT C{TQ ^ q)) (Sq), (145) 

where &c,(T Q (5q)) is the form a-q defined in the the formula (141) with q = C( t q(Sq))- The form a is a 
legitimate generating form since on fibres of the tangent fibration tq it is convex as a consequence of 
Proposition 6 and also positive homogeneous. 

Example 10. A point with configuration q' is tied to a fixed point go with a spring of spring constant 
kl . A second point with configuration q is tied to q' with a spring of spring constant k. The point q' is 
left free and the point q is subject to friction. As in Example 7, the internal configuration space is the 
affine space Q = M x M of internal configurations q = (q, q') modelled on 5 Q = V © V. The control 
configuration space is the space Q = M of controlled configurations q and SQ is the model space. The 
dual of SQ is denoted by F. The canonical projection 

n:Q^Q:q=(q,q')^q (146) 

is the relation between the two spaces. The virtual work function of the system is the form 



a : Q x S Q -> R : (q, q' , Sq, Sq') » k' \g(q' - go), Sq') + k(g(q - g'), Sq - Sq') + py/(g(6q),6q). (147) 
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Together with the projection r) it defines a generating family (a, rj) for the costitutive set. The critical 
set 



Cr(a, V ) = {q= («, </) € Q ; \j 5q , eV k'(g(q' - q ),Sq>) + k{g(q' - q), 5q') > o} 
= {<!= (<Z, q')eQ; k'g(q' - q ) + kg(q' - q) = 0} 



is the image of the section 



C : Q ->■ Q : q ^ (q, qo + - - - (q - qo j) 



of i]. It follows that the family (a, rj) can be reduced to the generating form 

k'k 



a:Q x SQ 
The constitutive set is the set 



.:(q,5q) 



k' + k 



(g(q - qo),Sq) + p\\s q \\. 



S ={(«,/) G Q x F ; ^^<<?(5 - «b) - /, Sq) + p\\Sq\\ > 



f)eQxF; 



k'k 



g{q - qo) - f 
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